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This  report  deals  with  one  aspect  of  restoration  of  atmospherically 
degraded  images;  the  problem  of  defining  isoplanatic  patch  size.  The 
general  restoration  situation  is  considered  to  indicate  the  situation 
where  the  concept  of  isopl anati cism  is  appropriate.  Then  a  model  is 
considered  and  some  problems  of  the  Fourier  transform  technique  of 
defining  isoplanatic  patch  are  pointed  out.  If  the  degraded  reference 
function  is  not  identical  with  the  point  response  function  then  iingo- 
lari ties  in  the  spectrum  of  the  intensity  of  the  resto^d  image  can 
give  rise  to  negative  image  functions.  A  measure  for  indexinq  this 
effect  is  indicated. 
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INTRODUCTION 


This  is  the  third  quarterly  technical  report  tor  Contract  No. 
f  30602- 72 -C-0305 .  The  object  of  the  program  is  to  provide  theoretical 
studies  of  selected  topics  associated  with  light-bean  propagation 
*h rough  a  turbulent  atmosphere.  The  present  problem  of  interest  has 
to  do  with  atmospheric  imaginq  and  image  restoration.  The  problem  is 
to  define  and  estimate  the  size  cf  the  ^soolanatic  paten  for  Atmospheric 
image  restoration  schemes.  This  problem  is  of  particular  importance 
to  those  intere_,r?d  in  taking  detailed  photographs  of  ear tn  sate " lues 
with  ground-based  telescooes.  It  mav  also  be  of  interest  to  loose 
doing  long  distance  imaging  of  terrestrial  objects. 

Imaqo  restoration  schemes  are  based  on  tne  concept  that  if  the 
degradation  for  a  single  point  were  known  it  would  give  sufficient 
information  to  reverse  tie  degradation  process  because  che  image  is 
basically  a  group  of  points,  (the  object  convolved  with  the  point 
rcsDonse  function),  and  it  is  assumed  that  all  points  are  degraded 
identically.  That  assumption  would  be  completely  valid  if  light  from 
every  point  on  the  object  and  from  the  reference  source  came  through 
exactly  the  same  atmospheric  index  fluctuations.  Since  the  object  is 
a  cluster  of  points,  the  question  is  then  how  far  can  the  points  on  tlte 
object  be  separated  from  the  reference  source  before  the  atmosphere 
affects  the  various  degraded  point  images  sufficiently  differently  that 
no  appreciable  image  restoration  can  be  accomplished.  The  linear  dimen¬ 
sion  of  such  an  area  at  the  ranye  of  the  object  is  the  isoplanatic  patch 
size.  It  provides  a  limit  on  the  size  of  the  object  that  can  be 
restored  with  a  single  rt  Cerence  source. 

This  report  discusses  some  of  the  schemes  for  image  restoration 
in  general  and  then  presents  preliminary  calculations  using  a  scheme 
capable  of  estimating  the  isoplanatic  patch  size  for  a  simple  situation. 
In  the  first  section  various  restoiation  schemes  are  indicated,  the 
concept  of  isoplanatic  patch  is  discussed  and  restrictions  to  the 
applicability  of  isoplanaticism  are  indicated.  In  Section  II  a 
particular  approach  is  introduced.  It  is  based  on  a  truncated  series 
representation  of  the  phase  of  the  incoming  electromagnetic  field. 
Several  other  possible  representations  are  also  ind  catcd.  Section  III 
contains  derivations  of  general  expressions  for  atmospherical ly 
degraded  images  and  for  images  restored  using  the  Fourier  transform  tech 
nique  These  expressions  are  then  discussed  in  the  light  of  the 
question  of  isoplanatic  patch  size,  and  of  certain  restoration  problems. 
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DISCUSSION  or  ISOPLANATIC  PATCH 


In  this  section  we  consider  situations  in  which  the  concept  of 
isoplanatic  patch  is  useful.  Various  points  arc  made  indicating 
restrictions  useful  in  its  definition.  Several  schemes  ltd ve  been 
suggested  for  the  restoration  of  atmospherically  degraded  images 
(Woods  Hole,  1 0f>6 ) .  Among  them  are  the  techniques  of  least  squares 
fitting  (Helstrom,  i%7),  ent^onv  maximization  (frieden,  1072), 

Fourier  transform  deconvuluti on  (Harris,  I  hth),  (I-C  iilann/>  l‘'*U7), 
holographic  phase  correction  (Gcoanan ,  1%6),  (fiast  ill,  'Hob).  .peckle 
interferometry  (Gezari ,  1072)  ar.d  point  bv  point  wave  front  pre¬ 
processing  (Vyce,  1973).  Injaae  processing  in  general  nas  also  been 
discussed  (Huang,  1971),  (Andrews,  1972).  Of  these  trie  Fourier 
trarsform,  holographic  and  preprocessing  schemes  are  of  interest 
because  they  obtain  their  information  from  an  additional  point 
reference  source. 

There  are  several  points  in  common  with  these  procedures  that 
help  to  define  the  isoplanatic  patch  study.  The  first  is  that  the 
Isoplanatic  patcli  basically  refers  to  instantaneous  single  recordings 
where  it  is  hoped  that  the  reference  source  contains  sufficient  in¬ 
formation  so  aj  to  allow  complete  image  restoration .  An  ensemble 
average  is  not  employed.  Indeed  if  the  ensemble  average  image  and 
reference  fields  were  to  lie  used  they  would  end  up  bcinq  independent 
of  position  because  of  assumed  atmospheric  homogenioty .  Short 
exposures  are  requi  red  because  temporal  averaging  effectively  brings 
into  action  too  many  ensemble  member;  with  similar  results, 

The  second  point,  although  well  known,  needs  documenting.  It 
concerns  the  telescope  aperture  size  requi rements .  If  a  telescope 
aperture  is  much  less  than  Fried's  coherence  length,  r0,  then  the  reso¬ 
lution  element  size  is  large  arid  motionless.  As  the  aperture  size  is 
increased  the  individual  image  elements  move  about,  but  maintain  shape. 
At  this  time  the  image  would  appear  distorted  so  that  a  grid  of  sharp 
straight  lines  would  appear  as  a  grid  of  sharp  but  wavy  lines.  In  such 
a  situation  image  detail  remains  but  relative  distances  must  be 
inferred.  As  the  telescope  aperture  becomes  still  larger  the  airy 
disc  pattern  becomes  deqraded,  i.e.,  broken  up.  It  seems  reasonable 
that  in  the  medium  aperture  situation  with  the  wavy  lines  that  a  fair 
amount  of  detail  is  still  available.  Further,  much  information  about 
relative  dimensions  can  be  inferred  by  trying  to  identify  lines  that 
were  obviously  straight,  such  as  construction  lines  or  edges. 

The  more  difficult  problem  occurs  when  the  point  response 
function  has  been  broken  up  so  that  detail  is  lost.  In  such  a  situation 
the  restoration  will  cause  the  most  gain  in  information  about  the  object 


If  the  distortion  of  straight  linr',  is  the  limiting  criterion 
then  possibly  the  isjplanatic  patch  could  be  estimated  from  angle  of 
arrival  correlations.  The  requirenent  is  basically  that  light  from 
the  reference  source  and  from  all  points  on  the  object  exhibit  the 
same  tilt  fluctuations  so  that  all  imane  points  move  simultaneously. 

The  typical  limiting  case  would  occur  when  light  from  two  point 
sources  had  arrival  angle  structure  functions  comparable  to  the 
diffraction  limited  uncertainty  indexed  bv  airy  disc  size.  This 
concept  has  not  yet  been  pursued. 

Another  point  concerns  the  do f ini  non  of  toe  isoplanatic  patch. 

It  should  be  closely  related  to  factors  describing  imuqe  quality. 

Two  of  these  have  already  been  sugpested  'Barry,  1971).  They  are  tne 
Rayleigh  resolution  of  a  restored  point  spread  function  and  the  Strehl 
ratio.  One  measures  the  resolution  of  the  restored  image,  and  the 
oJier,  using  the  relative  height  of  the  point  response  function,  gives 
an  indication  of  contrast  improvement.  Another  indicator  for  resolution 
might  be  the  percentage  decrease  in  integral  scale  width  under  the 
operation  of  restoration.  There  is  also  a  fidelity  criterion,  the 
normalized  mean  square  difference  uetween  object  ar.d  image  intensities. 

It  might  be  helpful  if  these  quantities  could  be  calculated  by  inte¬ 
grations  either  in  the  imaqe  plane  or  in  the  aperture  plane,  aperture 
plane  analysis  beinq  often  ccnvc-nie.it.  Thus  we  see  that  isoplanatic 
patch  applies  to  instantaneous  records  (photographs  or  holograms) 
where  a  point  reference  is  used  to  suoply  information  on  the  degradation. 
The  concept  is  most  useful  when  the  image  is  not  only  distorted, 
but  also  blurred,  finally  the  definition  of  isoplanatic  patch  should 
be  basically  couched  in  terms  of  how  the  image  is  degraded. 


II.  RLS1  ORATION  MOOLI. 

It  was  decided  to  investigate  the  restoration  process  based  on 
a  simple  model.  Since  the  process  is  in  essence  an  instantaneous  one 
an  approach  was  chosen  that  deals  with  the  instantaneous  field.  It 
was  also  decided  to  work  with  an  extremely  simple  situation,  that  is 
one  deqraded  imaqe  being  restored  with  a  degraded  point  reference 
source.  Further  it  was  decided  to  center  the  work  initially  on  the 
Fourier  transform  restoration  scheme.  The  last  two  items  are  in 
harmon)  with  experiments  already  in  progress  at  RADC. 

The  procedure  is  to  use  a  polynomial  expansion  that  could  be 
truncated  to  afford  maximum  simplification  Among  the  polynomial 
representations  considered  we;e  Fourier  series  representations  of 
either  the  aperture  plane  or  image  plane  fields  or  intensities 
along  with  the  cor^sponding  sampling  function  representation  of  the 
associated  transform  plane.  The  Gaussian-nermite  or  prolate  spher¬ 
oidal  function  representations  of  the  aperture  plane  or  image  plane 
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field  or  intensities  were  also  considered.  Further  a  Karhuner.-Loeve 
series  based  on  the  field  correlation  was  contemplated.  A  series 
expansion  of  tne  aperture  plane  phase  was  considered  and  chosen  because 
it  has  some  background  (Fried,  196b),  because  the  series  appears  to 
converge  after  a  very  small  number  of  terms,  and  because  the  series 
gives  a  chance  to  isolate  the  effects  of  wavefront  tilt  from  those 
of  higher  order  distortions.  This  represer  tation  assumes  negligible 
amplitude  fluctuations;  obviously  an  approximation.  However  it  is 
hoped  that  the  ensuing  simplicity  will  speed  the  results.  Fortner, 
the  situation  approaches  that  for  70.6  micron  radiation  at  tne  itADC 
propagation  range. 


Ill.  CALCULATION 

Wo  now  proceed  to  examine  on  a  somewhat  preliminary  basis  the 
process  of  restoration,  assuming  the  incominq  field  to  have  the  form 

(1  )  F  =  L'q  exp(j;). 

In  Eq.  (1) 


♦’>*.  Fn!x'y) 

where  the  polynomials  Fn(x,y)  are  equivalent  to  those  defined  by  Fried 
(Fried,  1965),  but  normalized  for  a  square  aperture  of  side  a. 


(2a) 
(2b) 
(2c) 
(2d) 
(2e) 
(2f ) 


F1  =  a" 

r  !jI  _x 

'z  "  "a  a/2 

F  ,  £j_. 
r3  a  a/2 


F  -  l3 

"  a  2  'I 


\  (x2  +  y2  -  (a2/6))/(a/2)2 


F5  =  1  (2/575/2 )(x2  -  y2 )/ (a/2 )2 


F6  =  (3/a)(xy)/(a/2)‘ 
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These  polynomials  are  chosen  orthonormal  so  that 


(3) 

/here 
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The  two-dimensional  spatial  spectrum  W(>  )  of  the  I  {■")  image- 
plane  intensity  is 


(4) 
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where  the  aperture  plane  coordinates  n  and  r2  have  been  replaced  by 
the  standard  sum  and  difference  variables 


and 

-  r2  +  rl 

n  “  - 2~~ 

and  7,  has  been  identified  with  ^ d/k  (d  is  the  image  plane  distance). 
In  order  to  evaluate  the  expression  for  W(<)  we  first  consider  the 

expression  E(o  -  E* ( n  +  Upon  inserting  the  expressions  for 

the  F^(r)  from  Eqs.  (2)  into  Eq.  (1)  and  simplifying,  we  obtain 
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where 

(6a ) 

A  - 
0 

=  2>3(d/a 

|?k)(a2.  x4a3.y) 

(6b) 

B 

XX 

*  6.-57? 

(a^  r  a5)  (d/2k  )2/  (a/2 ) 3 

(6c) 

B>.y 

=  Byx  " 

6a6(d/2k)2/(a/2)3 

(6d) 

Byy 

=  6/57? 

\  -  a5)(d/2k)2/(a/2)3. 

Jnserting  the  expression  in  Eq.  (5)  into  L'q.  (4)  and  performi  uj  the 
n  integration  gives 


(7) 


2  ,2  K  sinCKka/d)-  ,  '  XB,  .■  )] 

w(,.)  .  4E^(d/2,k)2  e  0 - k;B T.s 

xx  x  xy  y' 

sin[ ( (ka/d )- 1  )(B  -  +B v  ..  )] 


x  —  _ _ 

k  CB  »-  +B  <  1 
yy  y  xy  x' 


Equation  (7)  is  the  basic  expression  for  the  inctantareous  intensity 
spectrum. 

There  are  several  points  worth  noting.  First,  if  the  wavefront 
aberration  coefficients  (an  in  Eq.  (1))  go  to  zero,  then  the  expression 
reduces  to  the  triangular  function  as  expected.  Second  the  wavefront 
tilt  terms  a£  and  a3  act  only  to  shift  the  image  and  not  to  distort 
the  shape,  also  as  expected.  Third,  W(r)  does  not  factor  into  a 
function  of  ►  times  a  function  of  <».  This  arises  because  of  the 
term  BXy  which  is  proportional  to  the  distortion  coefficient  a6-  Indeed 
this  is  expected  because  Fs(r)  is  the  only  term  that  mixes  x  and  y. 


We  now  use  the  expression  for  the  degraded  image  spectrum  to 
give  an  indication  of  the  restored  image  spectrum.  Assume  that  two 
point  sources  are  imaged,  and  that  both  have  intensity  spectra  of  the 
form  in  Eq.  (7).  Then  the  restored  image  spectrum  W0(.  )  will  haye 
the  form 
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The  primes  indicate  that  tne  wavefront  distortion  coefficients  in  the 
Bjj’s  have  different  values  from  those  in  the  B{j's,  due  to  possibly 
different  atmospheric  paths. 

Nothing  has  been  said  about  the  definition  of  K0(- )  outside 
the  range  - (ka/d )• .  x,. y  +(ka/d).  It,  in  principle,  could  be  chosen 
by  analytically  continuing  the  two  pectra  W(.  )  and  W  (r)  and  then 
using  the  quotient  of  the  continued  parts  (Harris,  1971).  For  a  less 
laborious  approach,  it  could  arbitrarily  be  chosen  as  zero.  In  the 
situation  where  perfect  restoration  is  obtained  W0(«  )  =  1  for  x  , 

►  yj  ■  ak/d.  This  would  imply  a  (sinx)/x  function  for  the  restored 
image  intensity,  a  contradiction  since  the  intensity  can  not  be 
negative.  One  might  define  W0(T)  to  be  unity  everywhere  outside  the 
region  '*Xiiak/d,  .-y  ^ak/d,  giving  an  imoulse  function  for  the  perfectly 
restored  image.  However  for  imperfectly  restored  imaqes,  it  would  qive 
an  impulse  function  plus  a  second  function  no  matter  now  poor  the 
restoration,  which  seems  unreasonable.  Possibly  a  Gaussian  function 
exponential  or  a  Hanning  type  roundoff  in  both  the  x  and  y  directions 
would  provide  a  reasonable  first  compromise. 

The  function  W0(.  )  has  some  interesting  properties.  If  the 
distortion  coefficients  04- and  a^'-a^  are  respectively  equal, 
then  W0(r)  =  1  indicating  perfect  restoration.  If  do  not 

equal  the  respective  ajj'^a^  then  Wq(>)  will  have  singularities  where 
W'  (.-)  =  0.  Presumably  the  singularity  is  very  close  to  a  zero  if  the 
distortion  coefficients  do  not  differ  by  much.  Further  the  effect  of 
the  singularities  goes  to  zero  as  the  corresponding  zeros  approach.  How¬ 
ever  singularities  must  still  be  reckoned  with. 

The  effect  of  the  singularities  can  be  serii -qual i tati vely  con¬ 
sidered  by  examining  a  somewhat  similar  case.  Consider  the  function 
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then 


(1°)  F(W3)  =  ^  e  >x  W3(,-)d, 


J'0xO  +  ) .  -jx‘(l  +  -  )  _  •  o x 

e  n  dx'  e  r  sgn(x-x'  )o 

where  syn(x)  =  x/'xj.  The  "restored  image"  function  f  (W3)  in  Eq.  (10) 
has  one  desired  property,  namely  that  if  tne  two  functions  are  identical, 
so  that  the  roots  of  both  the  numerator  and  denominator  in  Eq.  (9)  are 
at  identical  locations,  then  '>  =  0  and  a  perfect  restoration  is  achieved. 
If  such  is  not  the  case,  the  "restored  imaqe"  is  imperfect:  it  has 
negative  intensity  neqions.  Tnese  are  bad  features  for  well-behaved 
intensity  functions.  Equation  00}  has  imaginary  parts  but  these  cancel 
for  the  function  W0(.)  because  tj  e  poles  a^  symmetrically  located. 

Although  the  results  of  Eq.  (10)  are  based  on  the  replaceinent 
function  in  Eq.  (9),  there  is  no  reason  to  believe  that  the  same  results 
will  not  be  obtained  using  the  restored  spectrum  in  Eq.  (8) 

Equation  (10)  describes  image  plane  behavior  and  is  therefore 
in  a  form  to  be  used  for  the  definition  of  isoplanatic  patch,  for 
example,  the  magnitude  of  the  second  term  in  Eq.  (10)  is  positive  and 
goes  to  zero  for  perfect  restorati on .  The  ratio  of  that  to  the  averaqe 
value  of  the  first  term  would  give  an  indication  of  the  relative 
effect  of  the  lack  of  isoplanatici ty .  The  coefficients  in  the  n  are 
proportional  to  the  wavefront  distortion  coefficients  (see  fns .  (6b)  and 
(0)),  and  might  be  estimated  by  using  the  mean-square  values  of  the 
distortion  coefficients  derived  using  turbulent  atmospheric  flu  tuations. 
"his  definition  will  not  be  recommended  particularly  at  this  tinv;  it 
is  only  used  as  an  example. 

The  discussion  presented,  although  preliminary,  does  present  an 
approach  to  the  problem  of  defining  isoplanatic  patch  in  the  presence 
of  some  problems  in  using  the  Fourier  transform  approach  to  restoration. 


IV.  SUMMARY 


This  report  has  presented  a  preliminary  discussion  of  some 
concepts  and  problems  associated  with  the  derivation  of  an  expression 
for  isoplanatic  patch  size.  First  the  situation  was  considered  in 
general  defining  the  case  where  isoplanatic  patch  might  be  a  useful 
concept  and  how  it  should  best  be  defined.  Then  a  particular  model 
for  an  atmospherically  distorted  wavefront  was  chosen  and  based  upon 
that  model  a  possible  approach  to  defining  isoplanatic  patch  was 
indicated. 
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